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Abstract A technique of timescale analysis performed directly in the time domain 
has been developed recently. We have applied the technique to study 
rapid variabilities of hard X-rays from neutron star and black hole bi- 
naries, 7-ray bursts and terrestrial 7-ray flashes. The results indicate 
that the time domain method of spectral analysis is a powerful tool in 
revealing the underlying physics in high-energy processes in objects. 
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1. Introduction 

The widely used Fourier analysis method in temporal analysis is to 
derive frequency spectra from a time series. After the Fourier transform 
of a light curve x(t)~) 



X(f j ) = Y,^k)e- i2 ^ kAt , 

k 

we can get the power density spectrum 

Pj = \X(fj)\ 2 



h = 3 it 



(1) 



(2) 



From two light curves, x\{t^) and ^(tfc), observed simultaneously in 
two energy bands at times and their Fourier transforms X\{fj) and 
X2(fj), we can construct the cross spectrum 



C(f j )=X* 1 (f j )X 2 (f j ) 



(3) 
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and then the time lag spectrum 

A(fj) = arg[C(/ J )]/2 7 r/ j (4) 
and the coherence coefficient spectrum 

r( f) = I < C ^ > I f5) 

KJ} V< \MfW><W^UW> 

People usually take a Fourier period 1// as a timescale and use the 
Fourier power spectrum to describe the distribution of variation am- 
plitude at different timescales, the time lag spectrum to describe the 
distribution of the emission time difference between two energy bands 
at different timescales and the coherence coefficient spectrum to describe 
the distribution of degree of linear correlation between high and low en- 
ergy processes at different timescales. But it is not correct to equate 
the Fourier period with the timescale; a Fourier component with a cer- 
tain frequency / of a light curve is not equal to the real process with 
the timescale 1//. For some important high-energy emission processes 
the Fourier spectra distort the timescale distribution of real physical 
processes seriously. It is needed to derive timescale spectra from ob- 
served light curves directly in the time domain without through Fourier 
transforms. 

2. Power Spectra 

The power density spectrum can be derived directly in the time do- 
main (Li 2001). For a counting series x(k), k = 1, ...,N obtained from 
a time history of observed photons with a time step At, the definition 
of variation power is 

Var(x) ^E£Li(x(fc)-x) 2 1 " 2 2 

P{At) = jAtr= (Atf = ^X>«-r) rms (6) 

where r = x/At. The power density p(At) at a timescale At can be 
then derived 

Fy ' &At At 2 -Ah ' K ' 

From (6), (7) we can calculate the power density for Poison noise 

T 

— X7 rms 2 /s (8) 

and the signal power density can be defined as 

Psignai(At) = p(At) - p noise (At) rms 2 /s (9) 
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Figure 1 Distribution of 
power density vs. time 
scale of a shot model. Top 
panel: the signal, stochas- 
tic exponential shots with 
time constant r between 
5 ms and 0.2 s. Mid- 
dle panel: simulated data 
which involved both signal 
and Poisson noises. Bot- 
tom panel: signal power 
densities . Solid line - 
power density distribution 
of time scale expected for 
the signal. Dashed line 
- excess Fourier spectrum 
from the simulated data. 
Plus signs - excess power 
densities calculated by the 
timing technique in the 
time domain for the simu- 
lated data. 

At (s) 

In Fig. 1 the bottom panel shows the expected power spectrum (solid 
line) for a signal consisting of stochastic shots (shown in the top panel) , 
the timescale spectrum of power density (crosses) and Fourier spectrum 
(dashed line) from the simulated data (shown in the middle panel) of 
the stochastic signal. As the rise and decay time constant of a shot is 
randomly taken from the range between 5 ms and 0.2 s, there should 
exist significant variation power in this timescale region. One can see 
from Fig. 1 that the Fourier spectrum significantly underestimates the 
power densities at the shorter timescales. 

Timescale spectra of power density can help to reveal the nature of 
physical process around compact stars (Li & Muraki 2002). Figure 2 
shows the timescale spectra and Fourier spectra of power density for a 
sample of accreting black holes (left side) and neutron stars (right side). 
The five black hole candidates demonstrate a significant power excess 
in the timescale spectra in comparison with their corresponding Fourier 
spectra at timescales shorter than ~ 0.1 s, but the two kinds of power 
spectrum in the studied neutron stars are generally consistent with each 
other. Assuming that there exist a stochastic process with a charac- 
teristic time of ~ 0.1 s in the black hole systems and that a significant 
variability of the accreting neutron stars comes from stochastic processes 
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with characteristic time constants much shorter than 
the observed power spectra for X-ray binaries. 



1 ms can explain 
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Figure 2. Power density vs. timescale of black hole binaries (left side) and neutron 
star binaries (right side). Crosses: Timescale spectra. Dots: Fourier spectra 



3. 



Time Lag Spectra 



Measuring the relative time delay between photons in two energy 
bands can give useful clue to understanding the emission mechanism and 
emitting region (e.g., Kazanas, Hua & Titarchuk 1997; Hua, Kazanas & 
Titarchuk 1997; Nowak et al. 1999). With the cross-correlation analysis 
one can calculate a time lag A from two light curves. Traditional cross- 
correlation method fails to calculate any time lag shorter than the time 
step At. In studying a complex process only a time lag A is not enough, 
we need to know time lags at different timescales, i.e. the timescale spec- 
trum A(At). Using Fourier analysis we can get a time lag spectrum. But 
as the Fourier technique is powerless for detecting the variation power 
of a stochastic process at high frequencies (short timescales) range, the 
Fourier analysis method is also powerless for detecting the time lags at 
high frequencies (short timescales). 
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The timescale spectrum of time lag can be derived with a modi- 
fied cross-correlation technique (Li, Feng & Chen 1999; Li 2001). For 
two light curves x\(k), x 2 {k) with a time step At the modified cross- 
correlation function at r is defined as 



where v(t) = x{t) — x, cr 2 (v) = J2 v2 > x (f) ls the counts in the time 
interval (t, t + At). Find the time lag A letting C(r = A; At) = max for 
different timescale At, we can derive the timescale spectrum of time lag, 



For comparing the two kinds of spectrum for time lags, timescale 
spectrum and Fourier spectrum, we construct two simulated time event 
series which consist of two kinds of random shot component : rapid 
components with typical timescale 1 ms and time lag 5 ms and slow 
component with timescale 0.01 s and time lag 0.02 s. Figure 3 shows 
that the timescale spectrum can represent the time delay distribution in 
the real process, but the Fourier analysis is powerless to detecting the 
time lags in the short timescale range. 



C(t; At) = ^2v 1 (kAt)v 2 {kAt + T)/a(vi)a(v 2 ) 



(10) 
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Figure 3 Time lags be- 
tween two series of events 
{t} and {t'}. Each se- 
ries consists of three com- 
ponents, {t} — {x} + 
{y} + {z},{t'} = {x'} + 
{y'} + {z'}. The first ran- 
dom shot component {a;} 
and {x 1 } have exponen- 
tially rising and decay time 
r = 0.01 s and x'(i) = 
x(i) + 0.02s. The sec- 
ond random shot compo- 
nent {y} and {y'} have ex- 
ponentially rising and de- 
cay time t = 1 ms and 
v'(i) — y(i) + 5ms. To- 
tal shot rate is 1000 cts/s. 
{z} and {z'} are two inde- 
pendent white noise series 
with rate of 100 cts/s. The 
series duration is 2000 s. 
Crosses are from the time 
domain technique, Circles 
from the Fourier analysis. 
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The timescale spectral method for time lag analysis is a powerful tool 
in revealing the characteristic of emission process. Applying the time 
domain technique, Qu et al. (2001) detected time lags of hard X-ray 
photons from Cyg X-l in the short timescale region (At < 0.01 s or 
/ > 100 Hz) with PCA/RXTE data for the first time, giving strong 
constraint on emission mechanism; Feng et al. (2002) revealed the time 
delays of soft photons in terrestrial gamma-ray flashes (TGFs) observed 
by BATSE/CGRO, strongly supporting the discharging mechanism of 
TGF production. Investigating temporal properties of GRBs by the 
timescale spectral method is in process. Figure 4 shows the spectrum 
of time lag of soft photons and the energy dependence of time lag for 
GRB910503 detected by BATSE. 



GRB 910503 
lime lags of soft photon (20-60 keV) 
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Figure 4- Soft time lags of GRB 910503. Left panel: Timescale spectra of time 
lag. Right panel: Time lag of 20-60 keV photons vs. energy of hard photons. 

4. Discussion 

The complex variability of high-energy emission shown in different 
time scales is a common character for X-ray binaries, super massive 
black holes and 7-ray bursts. The variability is caused by various phys- 
ical processes at different timescales. For understanding the emission 
process of high-energy photons it is necessary to know the variation char- 
acteristics in different timescales quantitatively, i.e. to derive timescale 
spectra from observed light curves. In addition to the methods of mak- 
ing spectral analysis for power density and time lag in the time domain 
introduced in last two sections, the algorithms to calculate timescale 
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spectra for coherence, hardness, variability duration, and correlation co- 
efficient between two characteristic quantities have been also proposed 
(Li 2001). 

There now exist two kinds of spectral analysis: frequency analysis 
and timescale analysis. As any observable physical process always oc- 
curs in the time domain, a frequency spectrum obtained by frequency 
analysis needs to be interpreted in the time domain. But a frequency 
analysis is based on a certain kind of time-frequency transformation. 
Different mathematically equivalent representations with different bases 
or functional coordinates in the frequency domain exist for certain time 
series, a Fourier spectrum with the trigonometric basis does not neces- 
sarily represent the true distribution of a physical process in the time 
domain. The rms variation vs. timescale of a time-varying process may 
differ substantially from its Fourier spectrum. Figure 1 shows that the 
Fourier transform distorts the power density distribution of a stochastic 
process at short timescales seriously. The timescale analysis performed 
directly in the time domain can derive real timescale distribution for 
quantities characterizing temporal property. In comparison with the fre- 
quency analysis, timescale spectra from the timescale analysis can more 
exactly represent timescale distributions of a real physical process, and 
more sensitively reveal temporal characteristics at short timescales for a 
stochastic process. Until now Fourier analyses of variabilities of X-ray 
binaries observed by various instruments are all failure to detect hard 
X-ray lags at the high frequency region of / > 100 Hz. Our simulation 
study shows it is an intrinsic weakness of the Fourier method, that more 
sensitive X-ray detectors of next generation can still not observe high 
frequency lags with Fourier technique. On the other hand, the timescale 
analysis can already derive time lag spectra at short timescales reliably 
from existing data. 
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